Complete stabilization and improvement of the characteristics of tunnel junctions by 

thermal annealing 
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We have observed that submicron sized Al-AlO^-Al tunnel junctions can be stabiHzed completely 
by annealing them in vacuum at temperatures between 350° C and 450° C. In addition, low temper- 
ature characterization of the samples after the annealing treatment showed a marked improvement 
of the tunneling characteristics due to disappearance of unwanted resonances in the current. Charg- 
ing energy, tunneling resistance, barrier thickness and height all increase after the treatment. The 
superconducting gap is not affected, but supercurrent is reduced in accordance with the increase of 
tunneling resistance. 
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Submicron sized metal-insulator-metal (MIM) tun- 
nel junctions have been used widely in applications 
such as radiation detectors single electron transistors 
(SET) [2], and tunnel junction coolers and thermometers 
[3| . They are also the central element in superconducting 
quantum bits 0, H, @, 0j 9 1 ^-nd magnetic junctions Q are 
developed for memory applications. However, there are 
often problems with stability of the junctions in room air 
conditions, the most noticable being aging, i.e. a slowly 
creeping increase of the tunneling resistance Rt with 
time. Instability of the junctions is an obvious drawback 
when considering practical and commercial applications 
of tunnel junctions. 

One possible reason for aging is the absorption of un- 
wanted molecules into the barrier. Direct knowledge of 
these processes is sketchy, although it is known that vac- 
uum conditions can slow down aging significantly. [To| 
Another possibility is the relaxation of the oxide barrier 
properties. Typically, tunnel barriers in MIM junctions 
are fabricated by room temperature thermal oxidation of 
a disordered base electrode. This procedure results in 
formation of a polycrystalline or amorphous oxide layer, 
in which atoms are not in their global equilibrium posi- 
tions. The relaxation from this "glassy" state could take 
months at room temperature. Moreover, in AlOa; barri- 
ers studied here, chemisorbed oxygen ions bound on the 
barrier surface during fabrication have also been shown 
to be a source of instabilities, [llj 

A typical remedy for finding a better energy minimum 
is thermal annealing, the artificial acceleration of the re- 
laxation processes by heating up the sample for some 
time and then letting it cool again. For technologically 
important Nb/Al-A10a;/Nb junctions, annealing treat- 
ments in gas atmospheres (He, N2, Ar, air) were shown 
to chan ge t he barrier properties, always increasing Rx • 
[13, [T^. Il4 [Tsj On the other hand, for magnetic tun- 
nel junctions, improvements in their performance were 
obtained with vacuum anneals. [l6l. [T7| 

In this letter we discuss how Al-A10a;-Al tunnel junc- 
tions can be completely stabilized, and their character- 
istics improved by thermal annealing in vacuum. The 
general trend of increase of Rt with the annealing is re- 



produced. The bias voltage dependence of the tunneling 
conductance reveal that the junction characteristics be- 
come more ideal after the anneal. In addition, the charg- 
ing energy, the effective barrier thickness and height all 
increase due to the annealing. Previous work on annealed 
Al-A10a;-Al junctions [l^ is partly inconsistent with our 
results, which is not surprising, as they used a forming 
gas (H2/N2) atmosphere in their process instead of the 
inert vacuum used here. 

A total number of 178 Al-A10a;-Al samples were an- 
nealed, containing either a single junction or two junc- 
tions in series, fabricated by conventional electron beam 
(e-beam) and standard double angle vacuum evaporation 
techniques. All the junctions had either a junction area 
~ 0.05/Ltm^ or ^ 0.15/im^. All samples were fabricated 
on oxidized silicon substrates, which were usually cleaned 
with an oxygen plasma just before Al deposition. For a 
few samples which were not annealed, plasma cleaning 
was not performed to see how the aging process is af- 
fected. The 300 nm wide and 55 nm thick Al wires were 
e-beam evaporated in high vacuum at a rate 0.3 nm/s, 
and the AlO^; barriers were thermally grown in situ in 
pure oxygen at ^ 0.2-0.5 mbar for 3 — 5 minutes. To re- 
duce the effects of contamination, all samples were post- 
oxidized in situ as a last step in ~ 0.5 bar of pure oxygen 
for five minutes. The samples were annealed in a vacuum 
chamber (P < 10^'^ mbar) to temperatures ranging from 
200°C to 500°C. The temperature cycle consisted of a 
quick (minutes) heating phase followed by a slow (hours) 
cooling cycle; details can be found in Ref. [Toi]. 

Aging was monitored before and after the annealing by 
measuring the sample resistance at room temperature vs. 
time, with all samples stored in room air. Representa- 
tive results are shown in Fig[TJ The solid and dashed 
lines show the typical aging of non-annealed samples 
with and without plasma cleaning of the substrate, re- 
spectively. The two samples shown were fabricated in the 
same evaporation and oxidation step to ensure identical 
conditions. It is clear that plasma cleaning of the sub- 
strate before Al deposition results in significant reduction 
of aging, but does not stabilize the junctions. We con- 
clude that some organic contaminants, most likely resist 
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FIG. 1: Aging of a non-annealed sample with (dashed) and 
without (solid) plasma cleaning of the substrate. Open circles 
represent aging when the sample was annealed at 200° C, open 
squares at 300° C and filled squares at 400° C, respectively. 
All resistances are normalized with initial resistance at time 
equals zero (7?o). 



residue on the substrate, are partly responsible for the 
aging phenomenon. The aging for a few representative 
annealed samples is also shown in Fig. [l] with different 
symbols denoting different annealing temperatures. The 
annealing process started approximately 20 minutes af- 
ter the end of fabrication, which defines the zero of time 
and the initial resistance Rq = R{t — 0) for each sample. 
Clearly, annealing temperatures 200° C and 300° C do not 
lead to stabilization, whereas complete stabilization oc- 
curs at 400°C. Note that the last datapoints correspond 
to times over two months after the fabrication. Complete 
stabilization was also observed at annealing temperatures 
350° C and 450° C (data not shown), but 500° C always 
resulted in breakdown of the junctions. In addition, the 
thermal treatments always led to increase of the sample 
resistance. At 400°C, where we have the largest amount 
of data (120 samples), the increase varied from a factor 
of two to a factor of eight, typically being between a fac- 
tor of two to four. The range of initial i?T /junction for 
the 400°C treatments was from 1 kfl to 100 kft, and no 
clear dependence of the Rt increase on the initial Rt 
was observed. Curiously, the low temperature anneals 
at 200° C and 300° C that did not lead to stabihzation 
show approximately the same logarithmic aging rate as 
the non-annealed samples, see Fig. [U but with a jump of 
the resistance caused by the treatment. This means that 
some of the resistance increase is not correlated with the 
stabilization in any way. 

For ten two junction samples, the normal state tun- 
neling properties were characterized before and after the 
annealing at 400° C by measuring the tunneling current 
and differential conductance dl/dV vs. bias voltage at 
4.2 K with standard lock-in techniques. In this case, the 
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FIG. 2: (Color online) A typical normalized dl/dV vs. V 
spectrum in units of Gt of a sample before (solid) and after 
(dashed, red) annealing, (a) shows a narrow bias range fo- 
cusing on Coulomb blockade, while (b) displays a wide bias 
voltage range (tunneling spectrum). 



low bias voltage regime (|T^| < 20mV) provides informa- 
tion on the charging energy, i.e. capacitance, of the junc- 
tions. The charging energy Ec = e^/(2Cs), where Cs is 
total capacitance of the junctions and the island, can be 
obtained from the height of the zero-bias Coulomb block- 
ade (CB) dip of dl/dV [Fig. U (a)] when fc^T » Ec, 
which is the appropriate limit for our junctions at 4.2 
K. For two junctions in series, the height of the dip AG 
and Ec are related by AG/Gt = Ec/i^keT) fs^, where 
Gt is the tunneling conductance around V — Q with- 
out the CB dip. Figure [2] (a) shows a typical CB dip in 
our samples before and after the annealing. Usually Ec 
increased by 5 - 10 % after the annealing, with initial 
Ecu ranging from 9.8 /ieV to 53^eV corresponding to 
capacitances Ce = 8.2-1.5 fF. The values of Rt — '^/Gt 
increased 1.5-3 times for these samples after the anneal- 
ing, varying initially from sample to sample between 9.5 
m and 151 kVl. 

In Fig. [2] (b) we plot a typical tunneling conductance 
spectrum (dl/dV) in a wide voltage bias range \V\ < 300 
mV. The spectrum shows many resonances (peaks) be- 
fore the annealing, indicating that initially our barri- 
ers are far from perfect. These resonances are usually 
caused by unwanted impurities within the barrier, al- 
lowing inelastic or resonant tunneling at some specific 
energies [l^. Almost all resonances have clearly dis- 
appeared after the annealing treatment, in addition the 
strength of the voltage dependence has been reduced. 
One possibility discussed before is that aluminum 
hydrates present in the disordered barriers dehydrate in 
the annealing, causing the resonances due to 0-H groups 
to disappear. In an idealized case, i.e. by assuming a 
trapezoidal barrier and that direct tunneling dominates, 
the voltage dependence of the tunneling conductance 
is commonly parametrized in the WKB approximation 
as G = Go[l -I- jV/Vo f], where = J^o/{e^md^) 
and Go = e^Ay/2m(j)o/{h'^d) exp{—2dy/2m(l>o/h). Here 
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A is the junction area, 0o is the average barrier height 
and d — {m* /mY^'^dox is the effective barrier thickness, 
where dox is t he p hysical thickness and m* the barrier ef- 
fective mass. [20] Within this model we conclude that the 
curvature parameter Vb increases in the annealing treat- 
ment. Together with the increase of Ec and Rt we can 
deduce that both the barrier height (/>o as well as its ef- 
fective thickness d increase with the annealing treatment. 
We stress that the increase of Ec alone (decrease of Ce) 
does not guarantee that d increases, as the dielectric con- 
stant most definitely also changes during the annealing. 
Similar kind of changes in junction parameters have also 
been reported earlier at least for gas annealed Nb/Al- 
AlO^-Nb junctions [13, The change in the absolute 
values of the barrier parameters are difficult to estimate, 
as the numerous resonances in the conductance spectrum 
before the annealing makes fitting to the model unreli- 
able. Quadratic fit to the annealed data yields the values 
d =9.5 Aand 4>q =1.6 eV for the sample in Fig. [21 

In addition to the improved tunneling characteristics, 
we also observed an improvement in the Al film quality 
by the treatments. This was checked by measuring the 
resistance of a 70 nm thick, 550 nm wide wire sample at 
4.2K before and after annealing at 400°C. The resistivity 
decreased by ~ 30% from 5.1 to 3.4 /zficm , corresponding 
to an increase in mean free path from 7.7 nm to 12 nm. 
High resolution SEM images of the film were also taken 
before and after the annealing, without any noticeable 
difference in the film grain structure. 

In addition to measurements at 4.2 K, we also deter- 
mined the superconducting properties of six unshunted 
single junction samples, cooled to below 100 mK in a 
dilution refrigerator. A representative I-V characteristic 
before and after annealing at 400° C is presented in Fig [3] 
(a). The increase of Rt is apparent as the decrease of the 
slope at high bias, in this sample approximately a factor 
of three. The superconducting gap A was determined 
from the curves, with the maxima of the dl/dV giving 
2A, whose value did not change after the annealing. Re- 
duction of the critical current Ic around zero bias is also 
visible, as expected from the Ambegaogar-Baratoff rela- 
tion, [2] where IcRt = const for a constant A =. To 
study the supercurrent branch more accurately, we mea- 
sured the switching currents by the cumulative histogram 
method, where a set of current pulses of constant height 
Ip were fed through the sample at a frequency 100 Hz 
and a duty cycle 5 %, and the time dependent voltage 
response was measured. By counting how often the junc- 
tion switched to the finite voltage state at a particular 
Ip, the switching probability distribution of the sample 
is directly determined. Fig. [3] (b) shows the obtained 
results for the sample whose IV data is shown in Fig. [3] 
(a). It is clear that the switching current Is (defined as 
the value where the probability = 0.5) is reduced to ap- 
proximately one fifth of the initial value, while the width 
of the distribution also decreases. 

By calculating the fluctuation free critical current Iq 



from the Ambegaogar-Baratoff relation and estimating 
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FIG. 3: (Color online) (a) Typical IV curve and (b) switching 
probability characteristics of a superconducting single junc- 
tion sample before (solid) and after (dashed, red) annealing 
measured at 80 mK. 

the capacitance of the junction as half of the capaci- 
tance measured for the two junction samples of the same 
junction size, the ratio of the Josephson coupling energy 
Ej and the charging energy Ec can be estimated to be 
Ej/Ec ^ 1.8. This means that the samples are in the 
underdamped regime, where for unshunted junctions the 
prediction is that Is oc Rrp^^^. ^1] The size of reduction 
in the switching current Is is therefore fully consistent 
with the increase in Rt- 

In conclusion, Al-AlO^^-Al tunnel junctions can be 
completely stabilized against aging by vacuum anneal- 
ing at 350 — 450°C . This is an obvious advantage, when 
considering practical and commercial applications of tun- 
nel junction devices. The improvement of the tunneling 
characteristics was also apparent as unwanted resonances 
disappeared after the treatment. Annealing led to a con- 
sistent reduction of junction capacitance and tunneling 
conductance for all samples studied, caused by an in- 
crease in the barrier height and thickness. Engineering 
these sample parameters in a controlled way by thermal 
vacuum annealing was not yet possible, partly because 
the microscopic understanding of the processes taking 
place is unclear. One possible application for the anneal- 
ing is to increase the sensitivity of Josephson junction 
threshold current detectors, currently used for example 
in superconducting quantum bit readouts or in shot- 
noise measurements. 22] With annealing, a reduction 
of the critical current can be achieved without a signif- 
icant reduction of the capacitance, which decreases the 
sensitivity. 23] 
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